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Abstract 



A Theory of Nuclear Magnetic Resonance (NMR) is developed for integer- 
spin, one-dimensional antiferromagnets, which exhibit the Haldane gap. We 
consider free boson, free fermion and non-linear cr-model approaches, all of 
which give similar results. Detailed anisotropy and magnetic field dependence 
is calculated and compared with experiment. 
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I. INTRODUCTION 



The one- dimensional Heisenberg antiferromagnet was argued to have a gap for integer 
(but not half-integer) spin by Haldanei, in 1983. Since that time a considerable amount 
of theoretical and experimental work has been done on this system. Many properties of 
this sytem have been calculated using relativistic quantum field theories in (1+1) space- 
time dimensions. While the non-linear a-model is probably the most firmly established 
model for this system, only a few of its features can be calculated exactly, so other more 
approximate models have been used extensively, notably free boson and free fermion models. 
All three models are generally expected to give qualitatively similar results. These models 
have been used to analyze neutron-scattering, susceptibility, electron spin resonance and 
other experiments (see [ and references therein). Our purpose here is to develop this 
theoretical framework to study nuclear magnetic resonance and to compare the theory to 
published experiments. Some of our results were obtained independently by Jolicoeur and 
GolinelliB. 

We consider the general Hamiltonian for S = 1: 

H = Y,{JS, ■ + DiS^Y + E [{S^f - [Sff] -ixbH-G- 8^}. (1.1) 

i 

For sufficiently small D/J and E/J, the model is still in the Haldane gap phase. In the pure 
Heisenberg model {D = E = 0) the lowest excited state is a triplet, at energy A ^ .4105J 
above the singlet groundstate. The D term splits the triplet into a singlet and doublet, and 
the E term further splits the doublet. A uniform magnetic field lowers the energy of one 
element of this triplet; it vanishes at a critical field. He. If the Hamiltonian is invariant 
under rotations about the field direction (eg. E = and the field is in the z direction), 
then there is quasi-long-range antiferromagnetic order above He- If the Hamiltonian does 
not have this invariance, there is true long-range order above He- The gyromagnetic tensor, 
G, is often assumed diagonal although we consider more general situations (see Appendix 

00 

The NMR rate, 1/Ti, is given by: 

/oo 
dte-'''Nt^bij^jb ^ {s\0,t),S^{0,0)} > 

0^ ^ ^M^ib ^ g-£;„/T (^^ ^ ^u.^/T\^ <n\S\Q)\m><m\S^{{])\n> 

i,j n,m 

2TT5{En- Em-OJN). (1-2) 

Here A^^ is the hyperfine coupling between the nucleus and the Ni spin, h is the direction 
of the RF field. In NENP, we will generally assume that A°'^{q) is not strongly g-dependent 
and is independent of h after averaging over the various nuclei. <> denotes a thermal average 
in an applied magnetic field, H (normally perpendicular to the RF field), cjtv = ^^NH is the 
nuclear Larmor frequency. This is of order ImK, negligibly small compared to the other 
relevant energy scales. Thus we are concerned with transitions between states n and m with 
essentially no change in energy but arbitrary change in wave-vector. We also can essentially 
set e^^l'^ = 1. We will concentrate on calculating 1/Ti in the low T limit, T <^ A, the 
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Haldane gap, where the low energy field theory approximation may be used. In this limit, 
there are contributions to 1/Ti from q k, tt and g ~ 0. 

There are several facets to a calculation of the NMR rate. We must consider contributions 
from two different ranges of q. Furthermore, the nature of the excitations changes radically 
near the critical field. As we vary the field at fixed temperature, we may identify three 
different regions: iib{,Hc - H) :§> T, fislH - Hc\ of 0(T), and fisiHc - H) <^ T. The 
behaviour of 1/Ti is quite different depending on whether or not there is rotational symmetry 
around the direction of the static field. Finally, the presence of defects, ie. chain ends, may 
lead to additional structure at temperatures well below A. 

In this paper, we mainly focus on the regime ^b{Hc — H) ^ T, A ^ T. In the rotation- 
ally invariant case, we give an asymptotically exact expression for the rate which is quite 
independent of any detailed assumptions underlying the field theories. Without rotational 
symmetry more complicated and somewhat model-dependent expressions are obtained. In 
both cases g ~ dominates at low T. Our result in the isotropic, zero field limit, was 
obtained previously by Jolicoeur and Golinelli. 

We also give a partial analysis of the region where /is(-ffc — H) <^T. We show that the 
free fermion model gives an exact description of this region, for q ^ 0, with or without axial 
symmetry. We also consider g ^ vr in this field region, showing that it actually dominates 
over q ^ 0. 

Where possible we compare our results to the experiments of Fujiwara et al.i on 
Ni(C2H8N2)2N02(C104) (NENP). While the overall agreement is fairly good, one puzzling 
discrepancy is encountered. 

In Sec. II we briefiy review the field theory treatment of this system and estimate the 
size of the q ^ n and g ~ contributions to 1/Ti. In particular, we review and generalize 
arguments that the free fermion model should become exact near the critical field. Sec. Ill 
considers the contributions to 1/Ti below the critical field for the isotropic and anisotropic 
cases. Sec. IV partially analyzes the region near criticality. Appendix |^ contains the details 
of the model calculations. Appendix ^ suggests a possible additional contribution to 1/Ti 
at low fields and low temperature. We expect the presence of non-magnetic defects, which 
essentially break up the system into finite length chains, to have a radical effect on the low 
energy properties of the system. Essentially decoupled 5 = 1/2 degrees of freedom appear 
at each chain end. This mechanism is not in itself sufficient to fully explain the low field 
experimental data. Finally, appendix makes note of structural properties of NENP which 
are expected to have effects on experiment. 



II. OVERVIEW 

At least three different field theories have been used extensively to discuss Haldane gap 
antiferromagnets. The first, and probably best founded is the non-linear u-mo delS. The 
antiferromagnetic order parameter, is represented by a three-component field, 0, obeying 
a constraint, 0^ = 1. The uniform magnetisation density is represented by the conserved 
density, /, of the field theory: 

l = {l/vg)(fxd(f/dt. (2.1) 
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Thus the spin operator at site i is written: 



Si^s{-iy(f) + l, (2.2) 
and the Heisenberg Hamiltonian density as: 

This is expected to give the asympotic low-energy long-distance behaviour of the Heisenberg 
model at large, integer-valued spin magnitude, s. The spin-wave velocity and coupling 
constant are given (up to higher order in 1/s) by: 

V = 2Js,g = 2/s. (2.4) 

Due to strong quantum fluctuations in one dimension, the groundstate is a (disordered) 
singlet and there is a gap, A to the lowest excited states, a triplet. For large s, 

A oc Je^^^ (2.5) 

(In the half-integer s case, an additional topological term must be added to the Lagrangian, 
leading to gapless behaviour.) The staggered spin field, 0, acting upon the groundstate, 
produces the triplet of single-particle states. However, the uniform magnetisation field, /, is 
a two-magnon operator, producing or annihilating a pair, or else fiipping the polarization 
of a magnon. Hence the low temperature neutron scattering cross-section is dominated by 
single magnon production for wave-vector q near vr and by two-magnon production for q 
near 0. 

Unfortunately, relatively little exact information is available about the cr-model. Ex- 
act S-matrix results indicate that the spectrum consists purely of the triplet of magnons, 
together with multi-particle scattering states but no boundstates. The exact matrix ele- 
ment of / between vacuum and two-magnon states is also known and can be used to predict 
the T = neutron-scattering cross-section at small q and energies between 2A and 4A 
(where 4-magnon production becomes possible). It is possible to add additional terms to 
the Hamiltonian, refiecting the crystal field anisotropy: 

n^n + a{<p'f + b[{<f)'^f -{<f)yf]. (2.6) 

However, no exact results are available on this model. Results are available from the large-n 
and related self-consistent approximations. In these approximations, the fields, 0", essen- 
tially become free fields with self-consistently determined masses and repulsive interactions 
which are assumed weak [of 0(l/n)]. 

The free, or weakly interacting boson model, sometimes referred to as the quantum 
Landau- Ginsburg modeli, is obtained by relaxing the non-linear constraint, 0^ = 1 and 
writing a simple 0^ Hamiltonian: 
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Here the are phenomenological mass terms, giving the gaps, which can be given different 
values for each of the three components, 0", to model crystal field anisotropy. A phenomeno- 
logical 101^ term is included, when needed for stability, but is treated perturbatively. 

A third, free (or weakly interacting) fermion model has quite a different motivation and 
is not obviously related to either of the above field theoriesB. The Hamiltonian is written 
in terms of a triplet of Majorana (ie. Hermitean) fermion fields, ipi and ipn- Here L and R 
label left and right-moving components. The Hamiltonian density is: 
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(2.8) 



Rather than being based on the large-s approximation, it appears to be special to the case 
s=l. It presumably becomes exact for the Heisenberg model with an additional isotropic 
biquadratic exchange interaction chosen to have the precise value to make the Haldane gap 
vanishS. (This model is on the phase boundary separating the Haldane and spontaneously 
dimerized phases.) This model has been shown to be equivalent to a = 3 Wess-Zumino- 
Witten (WZW) non-linear cr-model, or equivalently, three decoupled critical Ising models. 
The critical Ising model in turn, is essentially equivalent to a massless free Majorana (ie. 
hermitian) fermion. Reducing the biquadratic coupling moves the Ising models away from 
their critical point, corresponding to adding mass terms for the fermions. Four fermion 
interaction terms also appear, as written above, but these are generally ignored or treated 
perturbatively, in a similar phenomenological spirit to that of the Landau- Ginsburg boson 
model. In this representation, the uniform magnetisation density, / is again quadratic in the 
fundamental fields: 

1 = ^ [^L X V^L + ^R X ^r) . (2.9) 

The particles created by the fields, ip are identified as magnons with masses Aj. We see that / 
is again a two-magnon operator. It is also possible to represent the staggered magnetisation, 
in this approach, but it is considerably more complicated. Near the massless point, this 
operator reduces to the fundamental field of the WZW model, or equivalently to products 
of the order and disorder fields, a and pi of the three Ising models. These operators are 
non-local with respect to the fermion fields. The corresponding correlation functions can be 
expressed in terms of products of Pain-Leve functions. They exhibit poles at the fermion 
masses together with additional structure at higher energy. Unlike in the free boson model, 
a simple interpretation of the staggered magnetization density as a single magnon operator 
doesn't hold. 

To summarize, all three models have a triplet of massive magnons. In the first two models 
these are bosons and in the last fermions. In all models the uniform magnetisation density, 
/, is a two magnon operator. In the free boson or fermion models this is exact; including 
interactions in any of the models we expect some contribution from 4 or more magnons. In 
the free boson model the staggered magnetization, 0, is purely a single magnon operator. 
Including interactions there is an admixture of 3 and more magnons. Recent numerical work 
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indicates that this admixture is very smallO. This operator has a comphcated representation, 
non-local in the fermionic magnon fields, in the fermion model. 

We will be concerned with adding a magnetic field to the Hamiltonian. If the associated 
magnetization is conserved, then the effect is rather simple. This would occur, for example, 
if the field is in the 3 direction and rotational symmetry about the 3 axis is a good symmetry; 
ie. Ai = A2. Consistency demands that, in any of the field theories, we couple the field to 
the conserved density, /: 

n ^n- gf^sH -L (2.10) 
If the associated magnetization is not conserved, then extra terms could be addedi. We 



will use the term in Eq. ( |2.10| ) for simplicity. In the case where the magnetization is 
conserved, the energy of any state is simply shifted proportional to its magnetization, m, 
E E ~ gfiBHm. In particular, the magnon energy gaps become: 

MH) = A3 

A±iH) = A^±gfiBH. (2.11) 

Here the last two magnon eigenstates are created by the fields 0i ±202- If the magnetization 
is not conserved then the field dependence of the energies becomes model dependent. Ex- 
plicit formulas for magnon energies, Ea{k, H), have been derived in free boson and fermion 
mo The fermion result seems to agree very well with experiment; the boson one less 

so. As the field increases one of the gaps decreases, eventually reaching zero at a critical 
field of ~ 0{A/gfiB)- 

At low temperatures, we expect the nuclear magnetic relaxation to be dominated by q 
near or g near vr since this is where the lowest energy processes occur. It is clear that 1/Ti 
will vanish exponentially as T — since there is a gap. We wish to investigate whether 
the dominant contribution comes from g or g ^ tt. Let us first use the boson model to 
consider g ^ vr. From Eq. ([1.2D, we see that: 

1/Ti oc ^A^M^^J^e"^"/^ (l + e'^^/^) < n\^'\m >< m|0> > 

b,i ri,m 

6{En-E^-ujN). (2.12) 

Note that the expression is diagonal in spin operators. This is because we only consider 
Hamiltonians that have at least Z2 x Z2 x Z2 symmetry. Any correlators which contain cross 
terms of spin operators will be odd under some Z2 operation and must therefore vanish. 

We assume that the magnetic field is well below He. As we have mentioned above, the 
staggered magnetization field, is a single magnon operator in the non-interacting boson 
model. Thus it only has matrix elements between states whose energies differ by a magnon 
energy, Ea{k, H). In particular, there is no matrix element with energy difference (which 
is essentially zero). Including interactions in the boson model (as for instance given by the 
NLcT model) there will be a contribution from g vr at finite T. The simplest process is 
shown in Fig. |I|. It involves a 0^ interaction such as occurs in the quantum Landau- Ginsburg 
or NLcT models. The vertical line represents the field 0. The incoming line from the right 
represents a thermally excited magnon of non-zero momentum, k, and energy 2A. (For 
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simplicity, we consider the isotropic model, at zero field.) The two outgoing lines to the left 
represent magnons at rest. This diagram gives a non-zero matrix element proportional to 
A/A^. Note however, that since the energy of the initial and final state must be at least 2A, 
there will be a Boltzmann suppression factor of e~^^/^ in Eq. ( |1.2|) . Thus the contribution 
to 1/Ti will be proportional to: 

1/Ti oc A^e-^^/^. (2.13) 

Including anisotropy and a finite field, there will be various contributions of this type. It 
remains true that all these processes are suppressed by exp[— 2Aj-^jj-^(if)/T], where Ay^^y^{H) 
is the minimum gap. It is also possible to interpret this calculation as giving a finite width 
to the single magnon state at finite T. However, this cannot change the conclusion that there 
is a double exponential suppression factor, contrary to the proposed model in Ref. ( ^). 
Let us now consider contributions to 1/Ti from g ^ 0. These are given by: 

/oo 
rftA'^M^V'^'^* < {r(o,t),r(o,o)} > . (2.14) 

There is now, in general, a contribution even in the non-interacting boson or fermion model. 
This comes from the terms in which annihilate one magnon and create one magnon. In 
the presence of anisotropy, the three magnon branches are split, so we must distinguish 
inter-branch and intra-branch transitions (see Fig. |^). Note that both are possible since the 
wave-vector need not be conserved in the transition. However, the Boltzmann factor will be 
smaller for the inter-branch transition. In this case the optimum situation corresponds to 
a transition from a lowest branch magnon with a large kinetic energy to a middle branch 
magnon with zero kinetic energy. The Boltzmann factor is then exp[— Aj^-^^^^g/T]. Intra- 
branch transitions, in order to approximately conserve energy, will either have essentially 
zero change in wave-vector (forward scattering) or else the wave-vector will be essentially 
reversed (backscattering) . Clearly the dominant process will be the intrabranch transition 
for the lowest branch near the minimum energy k ^ 0. The Boltzmann factor is now 
exp[— Aj^jj^/T]. This contribution will be discussed in detail in the next section. 

We now wish to discuss the behaviour near and above the critical field. Clearly once 
we are sufficiently close to He that the lowest gap is of 0(T), the previous analysis breaks 
down. The q k, -k part may become comparable or even larger than the q {) part. We 
also need to consider possible infrared divergences and multi-magnon contributions. The 
nature of the critical point was established in Ref. ( §). If we assume exact rotational 
symmetry around the field direction (which is approximately true in NENP when the field 
is in the chain direction) then the phase transition is in the two-dimensional xy universality 
class, corresponding to Bose condensation of the lowest branch of bosons. There is quasi- 
long-range order for the staggered moment perpendicular to the field direction above He, 
and an associated gapless mode. Without axial symmetry, the phase transition is in the 
2-dimensional Ising universality class; there is true long-range-order above He-, with the 
correlation length and inverse gap only diverging right at H^. 

Rather remarkably, the free fermion model gives the exact critical behaviour of the 
uniform part of the spin operator, with or without axial symmetry. An argument for this, 
in the axially symmetry case, was given in Ref. ( ^) [See also Ref. ( |I^).] This follows from 
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the fact that the many body wave-function for a very dilute gas of repulsive bosons is just a 
free fermion wave-function multiplied by a sign function to correct the statistics. It is also 
fairly easy to see that the free fermion model becomes exact without axial symmetry. This 
follows from the well-known equivalence of the Ising model and a free Majorana fermion 
near the critical point. We emphasize that the free fermion model is certainly not exact 
far from He (although comparison with experiment indicates that it works well). Magnon 
interaction terms such as that in Eq. ( ^.8| ) must be included in general. However, they 
become irrelevant near He- This is quite easy to see in the Ising case. Two of the three 
branches of Majorana fermions remain gapped at the critical point and can be integrated out. 
Various interaction terms will be generated in the effective Hamiltonian for the remaining 
gapless fermion. However, there are no relevant interactions possible for a Majorana fermion 
since r vanishes by fermi statistics. Thus we can easily calculate the g ~ part 

of 1/Ti right through the critical point. This is done in Section IV. The q ^ ir part is more 
difficult since it requires correlation functions for the 2-dimensional Ising order and disorder 
variables away from the critical point on a cylinder of radius 1/T. A more sketchy discussion 
of this contribution is also given in Sec. IV. It turns out that the q ^ it part dominates very 
close to He- Well above He, in the axially non-symmetric case, when the gap becomes larger 
than T again, it is simpler to use the weakly interacting boson model, as discussed in Ref. 
( 1^), to calculate the staggered part. 



III. NMR RATES BELOW He 

A. T^^ With Axial (or SU{2)) Symmetry: Universality of Result, Exact Sigma 

Model Result. 

We first consider the isotropic case with the uniform magnetic field in the z-direction and 
the RF field in the x-direction. If the hyperfine couplings are also isotropic, A"-^ oc 5"'', then 
the intrabranch process doesn't occur. This follows because the RF field is perpendicular 
to the static magnetic field, so 1/Ti involves only matrix elements of S^. Since the three 
magnon branches have definite values of = —1,0,1, can only produce interbranch 
transitions. For more general hyperfine couplings, there will also be a contribution from S^. 
Defining, 

5,^ = l^e-''^^^(a;) (3.1) 

we see that in this case the leading low temperature relaxation rate will be determined by 
< k + q, —\Sq\k, — > and < k + q, —\Sg_^2k\ >• 1^; ~ > denotes the state on the lowest 

branch with wave- vector k and = 1. Here q must be very small since u!-{k + q)—u!-{±k) = 
ujn. At low T, k must also be quite small since the Boltzmann factor is exp(— tt;_(A;)/T). 
Since the total z component of spin remains a good quantum number in the presence of the 
magnetic field, it follows that for k ^ 0, 

<k,-\S'{0)\k,->=l (3.2) 

where we've chosen the following measure for the one particle states, 
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<k\k' >=2TT5{k-k'). (3.3) 

(This is true for arbitrary field, H.) In this case we can give a model independent, exact 
result for the low T relaxation rate. We also need the dispersion relation for the lowest 
branch at small k. By analyticity in k this must have the form: 

u_{k) = A-h+{vkf/2A + 0{k^), (3.4) 

where A is the zero field gap of the = ±1 branches and we may regard this formula as a 
definition of v; for convenience, we define 

h = gfisH. (3.5) 

In the Lorentz invariant field theory models, we have: 

uj.{k) = + {vky - h, (3.6) 



consistent with this form. Thus we obtain from ( p..2|) the following expression for intrabranch 
relaxation along the lowest branch: 

(V^Olntra « ^^lA^'l' J ^^(^-(^ + q) - ^-(k) - uj^) 

X e-"-('')/^| <k + q, -\S'{0)\k, - > p (3.7) 



TTf 2 



k^ + 2ujnA/v^ 



The square root in the denominator of Eq. ( |3.8| ) comes from integrating over the delta 
function. Note that it is necessary to keep at 7^ in this equation since otherwise there is 
a divergence at /c — 0, from the diverging density of states at the gap. 

At low fields, there will be interbranch transitions between the weakly split branches as 
well as intrabranch contributions from states with the "+" quantum number. The Boltz- 
mann suppression will always be that of the higher branch. The most important of these 
subdominant contributions comes from intrabranch transitions along the "+" branch. The 
reason is that low momentum transitions can occur in these processes causing a square root 
factor as in Eq. ( PTB| ) to be of order u^/v, ie. close to where the density of states diverges. 
The momentum transfer in the interbranch transitions must be essentially greater than the 
gap between the branches. The intrabranch transition rate along the "+" branch is identical 
to Eq. ( p. 81) , albeit with h —h in the exponential. 

In the regime of interest, cua? -C T <^ A, Eq. ( p.8|) gives: 



:Vri)intra « |A^1^^[log(4T/o.^) - 7](1 + e-2^/^)e-(^-'^)/^ (3.9) 



where 7 = 0.577216.... is Euler's constant. Essentially this formula, at h = 0, was given 
independently by Jolicoeur and GolinelliB, based on the large-N approximation to the NLa 
model. A similar expression was also derived by Troyer et. al. for the Heisenberg ladder 
problem having an excitation spectrum identical to the isotropic model discussed here0. 
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The present derivation shows that it is an exact result at T — > 0. We note that weak 
interchain couphngs, J', would also cut off this infrared divergence, replacing uj]^ in the 
above expression by a quantity of order V JJ'. 

Interbranch transitions will correspond to matrix elements 

I < A; + g, ±|^^(0)|fc, > |2 = I < A; + g, ±|^^(0)|A;, > 1^ ^ ^ (3.10) 

k and q are assumed small. The expression for the relaxation rate is identical to Eq. ( p.7|) 
with the appropriate matrix elements substituted and the energy of the higher branch in 
the exponential. Energy conservation enforced by the delta function will now replace a;7v by 
h{l + 2^) in Eq. ( p.8| ). One also needs to multiply the result by an overall factor of two 
corresponding to exchanging the labels of the states in the trace. The result is 

(l/TOMer « d^^l' + l^^1')^[log (^(ifx)) ' + (S-H) 

Writing: 

1/Ti = (l/TOintra + (im)lnter = F{h,T)e-^'^-^^/^ (3.12) 

following Fujiwara et. al.i; we see that F{h,T) decreases rapidly with h until h ^ T at 
which point the effects of the upper branches disappear. 

At somewhat higher T, the k dependence of the matrix elements will become important, 
as will the interbranch transitions. In the free boson model intrabranch matrix elements 
are, in fact, independent of h and identical for forward and backward scattering: 

< ±k, -|^"(0)|A;, - >^< ±k, -|/"(0)|A;, - >= 1 (3.13) 

The interbranch matrix elements (corresponding to transitions — ^ and + ^ 0) will be 

I < k,T\S^\k',0> P ^ - 
I )-ri I , I 2 



where ujk = \J (ykY + A^, u}'^, = \J {vk'Y + Ag and Aq is the gap to the 5*^ = branch (see 
the appendices for more details). Note that this result is still independent of the sign of k. 

In the free fermion and NLa models the expression for backward scattering is somewhat 
modified: 

< -A;,-|^"(0)|A;,- >= G(^) (3.15) 

The fermion model gives 

G{9) = -iA/uJk (3.16) 

while in the NLa model one gets 




(3.14) 
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(3.17) 



where the rapidity variable, 9 is given by 

cosh(^) = -1 - (3.18) 
Note that these are still independent of h. For vk <^ A, 9 = in + and 

G{9)^l-{{lr + \){vk/Ar (3.19) 

This is effectively due to magnon interactions and will, in turn, contribute a subleading T 
dependence to the prefactor, F{h,T). Estimating the strength of the magnon interactions 
from experimental measurements of 1/Ti would be difficult, due to various other sources of 
temperature dependence. 

Interbranch transitions in the fermion model are also sensitive to the sign of the scattering 
direction: 

2 

fti I < ±k,+\S-\k',0> 1^ (3.20) 

Notice that those interbranch transitions matrix elements which contribute to 1/Ti ^-re h 
dependent in all cases. This comes from the constraint uj±{k) = uJo{k'). 

Actually, all of the above arguments still hold with only U{1) symmetry when the field 
is along the symmetry direction. In principal, the function G{9) will remain independent of 
field but depends on how the SU{2) symmetry is broken down to C/(l). It is important to 
emphasize the universality of the leading term at at low T. The only requirement is that the 
system possess = ±1 degenerate low lying states with dispersion relation, cj ~ A + 
as A; — > 0. Once again, the addition of hS^ to the Hamiltonian will not change the matrix 
elements of P. The new energies will be a; ± /i, and with the low temperature assumption, 
one recovers the leading intrabranch contribution to 1/Ti. 



< ±k,-\S+\k',0> 



uj'u + k' 



+ 



I ujk±k 



B. ^ With Z2 Symmetry, Subcritical Magnetic Field Results 

Recall that for temperatures much smaller than the lowest gap, we need only consider 
the uniform part of the spin, /. We are thus led to calculate the transition matrix between 
two single particle states 

la,b{k:k')= r p-<k,afyt^O)\k',b> (3.21) 
J-00 Ztt 

where \k,a > and \k',b > are magnon states, and a and b denote the different magnon 
branches. 

The details of the calculation for the free fermion/boson models are given in the appen- 
dices. 
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Again, we perform this calculation in the low temperature limit {f3uj- ^1). In this 
regime, and with an anisotropic hyperfine tensor, we expect the dominant transitions to be 
intrabranch unless the branches cross or come close to doing so (this in fact happens when 
the field is placed parallel to the low branch axis). We start by calculating intrabranch 

contributions in both boson and fermion models. The relevant component of / is along 

the magnetic field. Taking this to be the 3 direction and the 1 direction corresponding to 
the RF field, Eq. ( |3.7| ) becomes 



Tf_i_ oc \A 



13 1 2 
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27r Jo 



(\ll_{k,k)\' + \ll^_{-k,k)\') 



uj^k) dk 
P + e{ky 



(3.22) 



The function e{k) is defined by 



(3.23) 



'dul 



k=Q 



and the three branches have the following dispersion relations 



(3.24) 



(A? + 



± 



4/l2 t;2^2 + 



(Af + Ai) 



(A?-Ai)2 



Bosons 



(3.25) 



(A? + Al 



± 



4/^2 t>2A;2 



:Ai+A2)2\ , (A2-A2)2 



U;2 = ^2^2 _^ ^2 



Fermions 



(3.26) 



When f3u!^ ^ 1 the integral in Eq. ( p.22| ) is strongly peaked at A; = 0. We can then 
make the following approximation (the validity of this will be discussed in a larger context 
later in this section): 



Tfi oc8|Ai3|2Ce-'^--(o)|/3_(o,o)|2i^f^(0)^ ' 

2n \ ok^ / 



(3.27) 



C= / dk , =- 
^0 + e(A;)2 



(3.28) 
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OO g 

dk— 



Using Eq. (|3.24| ) and making the change of variables P ^^r^ 2^71fi) ~^ ^' S^t 

C = I dk 



^ g^^{2.._(o)+^^)/3/4.._(o)^^ (c^jv(2cj-(0) + ujn)(3/Auj40)) (3.29) 

where i^o is the zero order modified Bessel Function. For u;_(0) ^ cuat, which is really 
the case we are considering, this reduces to 

C=-log(/5a;^/4)-7. (3.30) 



This is identical to the factor arising in Eq. (|3.9| ). 

For the boson model (see Appendix 



|f< (0.0)1^^4 '^■^^;>-V^^'''^--^^-> -^'^--^-'^/^ I (3,31) 

■ ' /2/i2(A? + Al) + (Af - Ai)V4 ' 



so F{h,T) goes to zero as a quadratic power of the field. 

lC(0.0)r^°|(|if)^ (3.32) 

Where we assumed without loss of generality, A2 < Ai. For the fermion model. 

This also vanishes with h"^. 

/I 7,2 A 2 

We can also extend the above approximations to intrabranch transitions along the "+" 
branch and interbranch processes. The former amounts to letting > +: 

oc 8\A're-'--''Y^A^M'^ (f^(0)) ' (- log(/5.;^/4) - 7) (3.35) 

The analogous expression for the interbranch rate is slightly more involved. Assum- 
ing that the branch labeled by r lies higher than that labeled by s, we can use the same 
approximations to arrive at 
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Q is defined by co'r(O) = uJsiQ). Notice that since S'^ is no longer a conserved operator, hence 
interbranch transitions + — are allowed. 
Once more we define 

Tf 1 = = F{h, T)e-^"-(°) (3.37) 

ab 

where one is careful not to double count in summing over a and h. 

There are regimes where it is no longer sufficient to calculate only one-particle transitions. 
When the temperature becomes comparable to the gap, multiparticle transitions become 
important. To include multiparticle transitions we simply replace the Boltzmann factor by 
appropriate occupation factors: /b(l + fh) = cosech^(^)/4 for bosons, and //(I — //) = 
sech^(^)/4 for fermions.0 At these temperatures it is also necessary to include the k- 
dependence of the integrand past the peak at the origin. We expect that at temperatures T a; 
y and fields /i ~ ^ the numerically integrated results would differ from Eqns. (|3.27| - p^36D 



by about 10 percent. Since we have not yet obtained exact expressions for the contribution 
of the staggered part of the spin to the relaxation rate and since, in the regime in question, 
single magnon contributions are expected to compete with those from the uniform part 
of the spin (this will be demonstrated in the next section), we will not bother to present 
numerically integrated results in this paper. 

In Figures |^-|^ we use Eqns. ( p.27| , p.35| and p.36|) to plot F{h, T) for bosons and fermions 



and for fields along the crystal a, h and c directions in NENP. We use A;, = 2.52mev, 
A(j = 1.17mev and Ac = 1.34mev, and assume a uniform hyperfine coupling to all spin 
components on a given site. We do this for varying temperatures and try to consistently 
account for contributions from relevant transitions including interbranch and multiparticle 
effects. Within approximations used, multiparticle effects amount to multiplying the final 
expressions by (1 ± e"^'^'*)"^. 

F{h, T) is shown for fields up to 9 Tesla even though the {(3uj- ^ 1) approximation is no 
longer valid at such fields. This is done to contrast the predictions of the boson and fermion 
models. It's easy to see that the boson result for F{h,T) diverges at the critical field, while 
no such catastrophe is present in the fermion result. This divergence is logarithmic and 
infrared. It will persist even after account is made for the staggered part of the correlation 
function. Multiparticle scattering will in fact worsen the effect. 

In NENP, when the field is along the b direction, we expect relevant interbranch tran- 
sitions only for small field. In this regime, one must also be careful to include intrabranch 
transitions in the second lowest branch. All these processes are of the same order. Even 
though the intrabranch rates vanish at low fields, the interbranch contributions are sup- 
pressed by the absence of low momentum transitions (ie. Q for the interbranch transitions 



is 0{y Af — A2) as opposed to O^ujn).) For this case, only need be calculated. 

When the field is along the c direction (corresponding to the middle gap), we restrict 
ourselves to calculating intrabranch transitions along the lower branch and interbranch ones 
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between the lower and c branch. There are no intrabranch processes along the c axis. 
Calculating the interbranch transitions amounts to calculating l/^cP ^^^1 |/acP- 

When the field is along the a direction, the calculation proceeds as above. The crossing of 
the branches provides for an interesting effect. At the subcritical field where branch crossing 
occurs there is an infrared divergence in the density of states (ie. e{k) oc k for small k). We 
expect this to be suppressed by the finite width of the applied uniform field or by a cutoff 
corresponding to interchain coupling but a peak should be seen in T^^ at this point. This 
peak can also be used to locate the true ac-axes for the chain (however, effects discussed in 
the appendices may broaden this peak considerably). 

As a brief aside, we would like to note that there are some curious differences in the 
behaviour of the boson matrix elements vs. that of the fermions. We have already seen that 
some transitions are fundamentally different in the two models even at low or zero fields. 
These may suggest experiments to further explore the nature of the low energy excitations 
in anisotropic Haldane gap materials. Even in the simple case of intrabranch backscattering 
the two theories differ (the fermion theory predicts a smaller, momentum dependent result). 
As another example, consider 

< 0,-|5j=o|0,3 > (3.38) 

ie. zero momentum transitions between the low branch and that parallel to the field (such 
a matrix element may be relevant in calculating ESR transition rates). The fermion model 
predicts a nearly field and gap independent result for such a rate in the limit /i ^ 0. 
The boson model prediction, however, is strongly dependent on the ratio of the gaps and 
therefore on the orientation of the field. To some degree, the reason that the theories feature 
such differences even in the isotropic case, is that the generator of rotational symmetry is 
fundamentally different in the two theories: the fermion magnetization does not couple 
left and right moving fermions, while the boson operator does. In particular, this means 
that intrabranch backscattering should vanish with the fermion mass, while no such miracle 
occurs with the bosons. 

There are some puzzling discrepancies between our result for F{h,T) and that of Ref. ( 
H). At about 5 Tesla, the experimental data gives Ft, < while our theoretical calculation 
gives the reverse at low temperatures. Moreover, the 'flatter' of the experimental curves 
is the one corresponding to the fleld perpendicular to the chain axis (about which there is 
approximate U{1) symmetry). The behaviour of the theories is quite easy to understand 
from the universal results valid in the axially symmetric case, discussed in Section (III. A). 
Fft is roughly independent of fleld with axial symmetry since is nearly h independent 
(in fact, -Ffo exhibits a logarithmic divergence as — 0). On the other hand, Fc vanishes 
quadratically as /i ^ 0. Including the small breaking of the axial symmetry corresponding 
to Ac — Aa = 2°K, Fh is essentially constant down to low flelds of order A^, — A^ ~ IT, before 
rapidly decreasing as seen in in the flgures. Given the original Hamiltonian, Eq. ( |1.1|) , we 
expect that our calculations are at least qualitatively correct. 

Finally, we would like to mention some recent NMR data collected on the 1-D S = 1 
spin chain AgVP2Se by Takigawa et. al.0. This material is highly one dimensional with 
a large gap (A ~ 320° K) and very nearly isotropic {6 ~ 4°K). These characteristics 
make it ideal for analysis using our results. There are, however, some questions about the 
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properties of the material which would have to be analyzed before an understanding of the 
NMR results is possible within the framework proposed here. The gap deduced from studies 
on the Vanadium atom (A ~ 410°i^) conflicts significantly with those performed on the 
Phosphorous sites and with neutron scattering data. In addition, the material has very 
low symmetry (corresponding to the space group Pljd) and very little is known about the 
possible small E and D terms in the Hamiltonian and their corresponding symmetry. There 
is fair qualitative agreement between the ^^P NMR data and our theory, and it is possible to 
explain some of the discrepancies using a temperature dependent anisotropic gap structure, 
but we feel that not enough is yet understood about gross features of the material to justify 
such speculation at this time. 



IV. CLOSE TO THE CRITICAL FIELD 

In this final section we would like to demonstrate that the staggered correlator con- 
tribution becomes crucial as h ^ he from below. First, we would like to describe the 
transition rate due to the uniform part of the spin-spin correlation close to the critical field, 
he — h <^ khT. In this limit, the long wavelength modes dominate the physics, and only 
intrabranch processes need be considered. Since the fermion model becomes exact in this 
limit, we will rely on its predictions. In the Z2 case, the dispersion relation becomes 

uik, h) = ^vlk^ + A2 (4.1) 

where vl = ^^^§7^^ and the gap is = {h - h,)^AiA2/{Ai + As)^ 

The U{1) case is different: the dispersion relation in the long wavelength limit is 

u;{k,h) = \ihe~h) + ^\ (4.2) 

In the fermion model, the coupling of the magnetic field to a conserved charge serves as 
a chemical potential for the magnon states. At criticality the dispersion is quadratic as 
opposed to linear as in the Z2 case. Above he the chemical potential drops below zero 
indicating that the ground state is now magnetic and filled with quasiparticles. 

In calculating the uniform spin contribution we must now include the appropriate occu- 
pation numbers in Eq. (|3.22| ): 



In the U{1) case we get a similar equation to (|3.^ ) but with appropriate occupation 
factors. The behaviour is logarithmic with temperature aX T ^ ujn even at the critical 
point. 

In the Z2 case we get 

{^^-% -f- ^7X^ roo^^(. + a..)cosech-(t) 

V 1 yUnif nv^ (Ai-A2)2io ^{u + unY-AI ^ ' 
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At criticality, we set u o: k and Ag = 0. We may simply rescale the integration variable 
to obtain 



(1/Ti)unif oc T (4.5) 

in the limit that Un T. This is expected from the Ising model where the uniform part 
of the spin corresponds to the Ising energy density operator, e, of scaling dimension 1. In 
terms of Majorana fermions this operator is ipL'4'R = ^• 

We can also say something about the behaviour of the staggered part of the correlation 
function at criticality. In both U{1) and Ising cases we know the critical behaviours of the 
staggered spin correlators. In both cases, the magnetic field acts as the temperature (and 
the temperature acts as the Euclidean time interval in an analogous 2D classical system). 
On the infinite Euclidean plane these correlators areB: 

< ^\z)(f){0) >~ l/\z\^2 (4.6) 

< a{z)a{0) >= l/\z\-i (4.7) 

The field = 0^ + is the charged U{1) field of the boson model; a is the disorder field 
of the Ising model (highly non-local in fermionic language) - it is the only local non-trivial 
primary operator aside from the energy density e. To get the finite temperature result, we 
simply make the conformal transformation from the plane into the finite strip0, Zp = e^'^*^"/^ 
to get 



< a{z)a{0) >-- 



sin(7rz//5)|' 



< 0t(^)0(o) >= .^""Z^]' (4.8) 



To get the contribution to Ti ^ we integrate over / dte^"^^^ and set z = it + e. 



Staff °^ / ^^ ^1, . , /^Mi (4-9) 
stag j_oo p^\sm{']Tz/P)\^ 



X is either 4 or 2. 

Since this is analytic as cjjv 0, we get that in the experimentally important limit 

Tf ^) g^^g cx + 0{ujn) Ising case 

%tag °^ ^'^^^ + ^^^^^ ^^^^ ^^^^ ^^-^^^ 

In the low T limit, this will be significantly stronger than the contribution from the 
uniform part of the correlator in both theories. In fact, as long as we are sufficiently close 
to the critical regime the above results will only be suppressed by factors of order 0(A/T), 
where A is the gap and A <^T. Thus in this regime, even slightly away from criticality, we 
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see that the most important contribution to 1/Ti comes from q = ti, ie. the staggered part 
of the spin correlation function. 

Farther still from criticality, the analysis breaks down but we expect the staggered spin 
contribution to influence 1/Ti through the region T = 0(A(/i)). 
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APPENDIX A: DIAGONALIZING THE BOSON/FERMIONS MODELS 

To do the necessary calculations we need to have a basis of eigenstates for the Hamiltonian 
and know the expansion of the field operators in terms of creation/annihilation operators for 
these states. Obtaining this is tedious (especially when the field does not lie in a direction 
of symmetry, for then all the branches mix) but the idea is to find the right Bogoliubov 
transformation which diagonalizes H. We do this first for the Boson model. 



As discussed in [ we seek to find the eigenvectors of the (non-hermitian) matrix rjM^ 
where 

Hk = d\Adk + aLfcA*a_fc + d\'Ba^_f^ + afcB*a_fc (Al) 



_T.*_A* (A2) 




We will now find the eigenvectors in the case where the field lies in a direction of sym- 
metry. In this case, only the excitations transverse to the direction of the applied field mix 
and we need only solve a (4 x 4) set of equations. Considering the two mixed components 
of 0, we write 

0Xfc,t = O) = -^[aLt^ + <] (A3) 



Ak.t = Q) = i^[a%-a^,] (A4) 
[4, ail] = 2Tr5ij5{k - k') (A5) 

ujq is arbitrary. 

Plugging this into the Hamiltonian, we then have that 

A='^l+^-\ha, (A6) 
4 Aujq 2 
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B = -^I + ^ (A7) 

where (T2 is the usual Pauh matrix, and 

K = dmg{Al + v^e,Al + v^e) (A8) 
Summarizing the above, we need to solve: 

(^OT _l_ K _ 1 L _ _'^0T I K \ 

which can be manipulated to give 

The eigenvalues of rjM are already known, as they are the solutions to the classical equations 
of motion and come in pairs ±uj± (see eqn. |3.25| ). Furthermore, we need only work to find 

one eigenvector of each pair because if ^ y ^ is a right-eigenvector of rjM with eigenvalue 

f Y* \ 

u, then ( 1 is a right-eigenvector of rjM with eigenvalue —u (X and Y are themselves 

two-component vectors). 

The bottom rows give the following set of equations 

= (A? + v^k:')x±,i - uJoUJ±^±,i + ihuj^i±,2 (All) 

= (A^ + v^k'')x±,2 - c^o^±e±,2 - ihuj^i±,i (A12) 

where it turns out to be convenient to work with % = X + y and = X — Y . The top rows 
can be manipulated to give 

= (/i^ _ uj%)x±,2 - ihujo^^^i + t^o^^±^±,2 (A13) 

= -(/i^ - ujI)x±,i - ihuJo^±^2 - (^o(^±^±,i (A14) 
These can then be worked to give 

(A^ + v^k^ + ujI- h^)x±,i = 2uJoUJ±Ui (A15) 

(A2' + v^k^ + iul- h^)x±,2 = 2cuocu±^±,2 (A16) 

{h^ + Al + - ujI)x±,i = -2ihuJo^±,2 (A17) 
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(A18) 



Note that if we fix the phase of Xi to be real then Yi must also be real as X2 and Y2 must 
be pure imaginary. The normalization condition now allows us to solve for the eigenvectors 
which form the columns of the transformation matrix between the old and diagonal basis of 
creation/annihilation operators. This normalization condition is slightly unusual since r)M 
is not hermitean: X^X — Y^Y = 1. In terms of x and ^ this is 



X±,lC±,l - X±,2^±,2 — 1 



(A19) 



The solution is 



X±,i 



1/2 



(A? + v'^k^ + ujI- h?){h? + + v'^k'^ - ul) ^ 



(A20) 



(A21) 



X±,2 



2ihujQ 



(A22) 



Ai + ^;^A;^ + a;^-/i^ 
2ujquj± 



-X±,2 



(A23) 



We can define the matrices x a-iid ^ with the rows labeled by the eigenvalues (+, — ) and 
the columns labeled by the original masses (1, 2). One easily verifies that in the hmit /i — > 
these become 



-It 



U!0 



i 







(A24) 



While in the limit A2 Ai, 



X 



VA2 + -f;2A;2 



<^0 



1 -i 



2 VA2 + -t;2A;2 I 1 i 



(A25) 



It turns out that 



(A26) 



(A27) 
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Where the 6's are the operators which diagonahze H. 

So now we have all the equipment necessary to calculate la,b{k, k'). 



(A28) 



where we define the cross product matrix with the Levi-Civita symbol by S* = e^^''. The 
matrices, x ^5 ^ire now three by three with the inclusion of the trivially diagonal unmixed 



component (ie. the three component): X33 



^0 



and 



33 



y/Aj+v'^k'^ V '-'0 

We would now like to repeat this diagonalization procedure for the Fermion model. The 
free Hamiltonian is 



n{x) = 



(A29) 



Setting f = 1, we can write 



R 



°° dk 
2^ 



(A30) 



Jo Ztt ^ 



(A31) 



{alai^} = 2TT6,,6{k-k') 
the Hamiltonian density in A;-space becomes 



where 



M 



Ik — ihx iA 

—iA —Ik — ihx 



(A32) 



(A33) 



(A34) 




5. = T (A35) 



The idea now is to diagonalize this matrix and find the eigenvalues and eigenvectors. In 
other words, find the unitary transformation which diagonalizes H. We assume the field is 
in a direction of symmetry so that we need only diagonalize a 4 x 4 matrix. Given that 
the field is in the 3 direction, the eigenvalues for the mixed states are ±u± from eqn. p.26| . 
It may be more illuminating to write out M in a basis that is more natural to the U{1) 
problem. Using 



^R,k 
(^R,k 




(A36) 
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In this basis, M becomes 

Where A = ^1+^^ and 5 — The equations for the components of the eigenvectors 

possess the symmetries 

ui <->■ U2, U3 ^ U4,h ■(-^ —h (A39) 



ui Us, U2 M4, a; —uj (A40) 
where cu is the eigenvalue. After some algebra, 

2i5(k — cu) , , 

2iA(fc + a;) , , ^ , 

Using the normalization condition, 

Ekl' = l (A44) 

we set the phase of Ui to be real for positive eigenvalues; the above symmetries allow us the 
freedom to choose a convenient phase for ui corresponding to negative eigenvalues. 

Ml = 26{k + uj){uj^ + A^ - (A; + h)^ - 5^) ^ (A45) 
[46\k + oof {{00^ + A^-{k + hf - 5^f + 4A2(A; - oof) + 
((^ + hf -e + S'- A2)2((^2 + - + hf - d'f + 4A\k + uf)]'^ 

We define the 6x6 diagonalizing matrix with columns given by the eigenvectors as 

- , ul_ , , mL^^ , , mL^3 , ) ( A46) 

al = UXIP^ (A47) 
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\/2 



/II \ 

i -i Q 

V2 

1 1 

-i i 

Voo oov^y 



V 

V(Ti 



(A48) 



The diagonal operators, /^"^ are defined as: 



Cfc 

4 



(A49) 



Our freedom in choosing the phases for the eigenvectors corresponding to negative eigen- 
values allow us to write 



X 



R T 
T R 



(A50) 



The (i's and c's correspond to left and right moving fermions, respectively. This becomes 
clear in the limit Ai — > A2 — > 0. 

Some limiting forms of R and T are: 



R{h ^ 0) = - 



V 



0)2 +fc 
UJ2 








V '^3 



(A51) 



ijJ2—k 
L02 



0)2— k 








V 



u>3—k 

L03 



(A52) 



R{6 0) 



1 
72 







a;A+fc 

OJA 















0)3 



(A53) 



T((5 ^ 0) 



x/2 



OJA 







V 








OJA 








V 



u)3—k 

U)3 



The expression for the uniform part of the spin in terms of fermions is 



(A54) 



(A55) 
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We can write la,b{k, k') as 

la,b{k,k') = ~i J dp dq < k,a\{al + dp) 'ij {dq + d*)\k',b > (A56) 



where again E is the cross product matrix. 
After some algebra, we get 



iaA^. k)- i\ j^TyT^yj^ ^ R'^ a^V^^V a,T T^V^^V*T* + R^a,V^^V*a,R* ' ^ ' 



Each index of this matrix runs over six states; the first and last three correspond to right 
and left movers respectively. In the case of U{1) symmetry each set would correspond to 
states of definite spin. 



APPENDIX B: IMPURITY EFFECTS ON 1/Ti 

The NMR experiments done on NENP have exhibited curious behaviour at small fields 
(0-4 T) and low temperatures. In this regime, our calculations have shown that 1/Ti 
should decrease rapidly with decreasing applied field. Experiments, however, indicate a 
non-zero lifetime at h = that undergoes a minimum at about 2 T before showing the 
expected exponential behaviour at around 4 T. Moreover, the zero field lifetime increases 
with temperature.! 

If one attempted to fit the mid field data (which is expected to be in good agreement 
with our calculations) to our results and assumed that there is no other mechanism for 
relaxation, then the observed values would deviate from the calculated ones at low fields 
[h = .2 — 3 Tesla) by a factor as big as 300 for T = 1.4 Kelvin. 

Clearly the low field bump is an impurity effect. Dynamics due to dimensional crossover 
are ruled out by the fact that 1/Ti increases with temperature at the low fields discussed. 
A possible, yet not fully convincing contribution can come from the dynamics of the spin-i 
degrees of freedom at the ends of large but finite chains in NENP. This is discussed further 
in ref. [ ^. The essential physics is that the resonance energy of a spin-i at a chain end 
in a uniform field crossed with an RF field is broadened by the interaction with thermally 
excited magnons even at temperatures well below the gap. At the Larmour frequency of 
a Hydrogen nucleus and low magnetic field, the Fourier transform of the local spin-spin 
correlator (which is proportional to 1/Ti) is of the form:i 

oc e-^/^ (Bl) 

This is true in both the limits of long and short chains. In both instances this is not consistent 
with the weaker temperature dependence seen at constant low magnetic field. We conclude 
that the experimental results are due to a mechanism not intrinsic to the Haldane system 
or coupled to it. 
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APPENDIX C: THE STRUCTURE OF NENP AND EXPERIMENTAL 

RAMIFICATIONS 



Each chain in NENP is comprised of Ethylenediamine-Nickel chelates separated by nitrite 
groups. It is important to reahze, however, that two neighbouring Ni^^ are not equivalent; 
rather, one is related to the other by a tt rotation about the h axis. Also, the angle along the 
— Ni — O bond is not exactly vr, meaning that the Ni site is not truly centrosymmetric. 
Most importantly, the local symmetry axes of each Ni ion are rotated with respect to the ahc 
(crystallographic) axes. To demonstrate this we now note the coordinates of the Nitrogen 
atoms in the Ethylenediamine chelate surrounding the Nickel (placing the Nickel at the 
origin) 



Atom 


a (A) 


h{A) 


c(A) 


iV(l) 


2.053 (3) 


.162 (3) 


.338 (3) 


N{2) 


.619 (3) 


-.184 (3) 


-1.971 (3) 



The other Nitrogen atoms in the chelate can be obtained by reflection through the Nickel. 
One easily sees that projecting this structure onto the h plane yields symmetry axes (in the 
h plane) rotated ~ 60° from the ac-axes. This is shown in Fig. |^. The inclination of the 
local Nickel axes from the ahc system can be obtained by taking the cross product of the 
two Nitrogen vectors (ie. the normal to the plane described by the four Nitrogen atoms in 
the chelate: 

n= (-.06 (1),. 98 (1), -.11 (1)) (CI) 

From the above, it is clear that the local Ni 6'-axis makes a ~ 10° angle with the 6-axis, 
while the azimuthal angle in the ac plane is ~ —28° from c. The 10° tilt is roughly about 
the a' direction of the local symmetry axes. 

One may worry that the NO2 radicals may distort the local symmetry axes, but remark- 
ably enough, when projected onto the ac plane, the three atoms in the molecule sit on the 
c' axis. This reinforces our suspicion that the local symmetry axes are indeed what we have 
above. 

Next, we consider the whole space group of NENP. The most recent attempt to solve for 
the crystal symmetries has concluded that the true space group of the material is Pn2i 
this is a non-centrosymmetric space group with a screw 2i symmetry about the h axis, 
diagonal glide plane reflection symmetry along the a axis, and an axial glide plane reflection 
symmetry along c. Experimentally, attempts to solve the structure in Pn2ia have not been 
successful; rather, it seems that Pnma gives a better fit. The main difference between the two 
is the presence in Pnma of a mirror plane parallel to h at centers of symmetry at various 
locations in the unit cell, and two-fold screw axes separating these centers of symmetry. 
The reason for the experimental discrepancy is attributed to disorder in the orientation of 
the nitrite group, the perchlorate anions, and the existence of a local or pseudo center of 
symmetry lying very close to the Ni (thousandths of an Angstrom)^!. A crucial point is 
that both space groups share the axial glide planes along a, the diagonal glide planes along 
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C and the 2i screw symmetry about b. These generate a total of 4 Ni sites per primitive 
cell and two chains through each cell. The two chains are such that the Ni chelates on one 
are the mirror image of the other. Figure ^ shows a projection of this picture onto the ac 
plane. 

The presence of the 2i screw symmetry about each chain axis introduces staggered con- 
tributions to the local anisotropy and gyromagnetic tensors. This is because, as motivated 
above, these are not diagonal in the crystallographic coordinate system. The resulting spin 
Hamiltonian is 

H = JY,[S,- S,+, + ■ BS, - fiBSi ■ GH + {-iy{S, ■ dSi - fisS, ■ gH)] (C2) 

i 

We make the assumption that the symmetry of the mass and (7-tensors is the same 
(ie. that at each site they can be simultaneously diagonalized) . We can get the required 
parametrization for the ^f-tensors from high temperature uniform susceptibility datalll. This 
is based on the idea that at high temperatures the Ni atoms will behave as an ensemble 
of uncoupled spins (s = 1) with the same gyromagnetic tensor as in the antiferromagnetic 
case. With this in mind we get 

( Gc' cos\e) + Gb' sm\e) \ 

G = Ga' (C3) 

V Gb'Cos\e) +G,,sm^i9) J 

/ sm{e) cos{e) {Gb' - Gc') \ 

g= (C4) 

\sm{9)cos{9){Gb' -G,,) 00/ 

Here 6 ~ 10°, and Ga' = 2.24, Gb' = 2.15, G^' = 2.20 are the values for the local G- 
tensor that give the observed high temperature (/-tensor when averaged over the unit cell. 
Correspondingly, the mass tensors must have the following form: 

D^, 

D = I I (C5) 

Db' 



'-^(Db'-D,, 

I (C6) 

'-^{Db'-D,,) 

The parameters Da', Dc', Db' are to be fitted by experiment to the model used to describe 
the system. 
Using 

S^ = (-l)> + 0xn (C7) 
The boson Hamiltonian can now be written 
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H = I cia;[|n2 + |(|^)' + ^0-D0-/iBi? -0(0x11)+ (C8) 

^0 ■ d(0 X n) - ■ gH + A(0')'] 

The term containing d breaks the Z2 symmetry along the a' (lowest mass) direction. It 
will also renormalize the masses. The second effect can be ignored in the approximation 
that the 0^ term is ignored if we assume the masses are physical. The first effect leads 
to the presence of a static staggered field even below a critical magnetic field. A gap will 
always persist. The staggered field term will break the Z2 symmetry along the c' or h axis, 
depending on whether the field is applied in the b or c' direction, respectively. A static 
staggered moment will likewise appear due to this term. 

We would now like to discuss the effect of having two inequivalent chains per unit cell, 
with local axes different from the crystalographic axes. We label the two chains found in a 
unit cell of NENP 'chain 1' and 'chain 2' corresponding to the chains in the upper right and 
lower left corners of Figure |] respectively. The dispersion branches of chain 1 are given by 
(20) of [ 1^ only the field is a — 30° from the c' axis where a is the angle of the field from 
the crystallographic c-axis. Similarly, the dispersion branches of chain 2 are calculated with 
the field a — 150° from the c' axis. We can now graph these and compare with experiment. 

Experiments which average over signals, like susceptibility or NMR measurements, 
must consider their results an average of two different experiments (corresponding to the 
two different chains with their relatively different field configurations), on the other hand, 
experiments such as ESR, should show a separate signal for each chain. The NMR relaxation 
calculations performed in the main part of the text ignore these effects. This ought not make 
a difference to the qualitative conclusions. 

Figure ^ shows the dispersions for chains 1 and 2 (solid and dashed lines, respectively) 
when the field is 7r/3 from the crystallographic c axis in the ac-plane. This is an example of 
how transitions at two field strengths ought to be possible in the ESR experiment. 



Figure |T0| shows the resonance field versus orientation of field in the crystallographic ac- 
plane. The lower branch denotes transitions in chain 1 while the upper branch corresponds 
to transitions in chain 2. The transitions were calculated at .19 meV. corresponding to 47 
GHz In addition the experimental results of Date and Kindoisl are represented by the crosses. 
One immediately sees that the data does not compare well with the predictions based on 
the free model used above, for instead of following one of the branches, the experimental 
results lie between them. Furthermore, it seems unlikely that perturbations will cause such 
a significant shift in the resonance field. One sees that the discrepancy is ~ ±1 Tesla. One 



possible explanation is that since the ESR signal in [ ^ was also ~ ±1 Tesla in width and 
symmetric (in conflict with the predictions of [ |^)], the signal from the resonances in both 
chains was somehow smeared and interpreted as one single peak. Seen that way the model 
predictions are in good agreement except for the large field regime. One also has to keep in 
mind that the high-field boson dispersions are not accurate and therefore the predictions at 
larger angles could easily be .5 Tesla (or more) off the mark. 

To end this discussion on, we'd like to elaborate on a previously made statement regarding 
the assignment of masses to the local Ni symmetry axes. It's easy to see that switching the 
masses around is tantamount to a 7r/2 shift in Figure ^ (the fact that the gyromagnetic 
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constants are not the same in orthogonal directions will not change the ESR resonance 
graph much since the ratio of the gyromagnetic constants is 0.98). Redrawing Figure 10 
with this geometry misses the experimental results by 4 Tesla at and 90 degrees, where the 
two chain resonances coincide. This determines the proper labeling of the local symmetry 
axes. 
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FIGURES 
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FIG. 1. First non- vanishing contribution to relaxation due to the staggered part of the spin 
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Magnon Dispersions For Lowest Two Branches 




FIG. 2. Inter- vs Intrabranch transitions 
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Boson F(H,T): 

Contributions From Tine Uniform Part of the Spin 
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Fermion F(H,T): 

Contributions From The Uniform Part of the Spin 
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FIG. 7. Local and crystallographic axes projected onto the ac-plane in NENP 
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FIG. 8. A projection of the NENP unit cell onto the ac-plane showing two chains per unit cell 
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FIG. 9. Dispersions for the two chain conformations and sample resonant transitions for 
uniform field placed 60° from the crystallographic c-axis. 
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ESR Resonance in NENP 




